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ABSTRACT 



T Pyx is a luminous recurrent nova that accretes at a much higher rate than is 
expected for its photometricahy determined orbital period of about 1.8 hours. We here 
provide the first spectroscopic confirmation of the orbital period, P = 1.8295 hours 
(/ = 13. 118368 ±1.1 X 10~^ c/d), based on time-resolved optical spectroscopy obtained 
at the VLT and the Magellan telescopes. We also derive an upper limit of the velocity 
semi-amplitude of the white dwarf, Kl = 17.9 ± 1.6kms~^, and estimate a mass ratio 
ofq = 0.20±0.03. If the mass of the donor star is estimated using the period-density re- 
lation and theoretical main-sequence mass-radius relation for a slightly inflated donor 
star, we find M2 = 0.14 ±0.03 M©. This implies a mass of the primary white dwarf of 
Ml = 0.7±0.2 M0. If the white-dwarf mass is > 1 M©, as classical nova models imply, 
the donor mass must be even higher. We therefore rule out the possibility that T Pyx 
has evolved beyond the period minimum for cataclysmic variables. We find that the 
system inclination is constrained to be i ~ 10 degrees, confirming the expectation that 
T Pyx is a low-inclination system. We also discuss some of the evolutionary implica- 
tions of the emerging physical picture of T Pyx. In particular, we show that epochs of 
enhanced mass transfer (like the present) may accelerate or even dominate the overall 
evolution of the system, even if they are relatively short-lived. We also point out that 
such phases may be relevant to the evolution of cataclysmic variables more generally. 

Key words: stars: individual: T Pyx — novae, cataclysmic variables. 



1 INTRODUCTION 

Binary systems containing a Roche-lobe filling main- 
sequence star that loses mass to a primary white dwarf (WD) 
are r eferred to as cataclysmic variables (CVs, see .Warner. 
I1995I for an overview). T Pyx is a luminous CV in which 
the donor star is transferring mass at a very high rate onto 
a high- mass white dwarf, resulting in unusually frequent 
thermo- nuclear runaways (TNRs) on the surface of the pri- 
mary. Between the years 1890 — 1967, T Pyx has undergone 
five such nova eruptions, with an recurrence time of about 20 
years between the eruptions, and was therefore classified as 
a member of the recurrent nova (RN) subclass. However, the 
last eruption was in 1966, which means that T Pyx has now 
passed its mean recurrence time by more than 20 years. The 
eruptive behavior of RNe in comparison with classical novae 
is thought to be due to a high mass-transfer r ate in combi- 
natio n with a massive primary white dwarf (JYaron et al.l 
I2OO5I ). 
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[Patterson et al.l (| 19981 ) carried out an extensive pho- 
tometric study of T Pyx and found a stable, periodic sig- 
nal at P = 1.83 hours that was interpreted as the likely 
orbital period. This would place T Pyx below the CV pe- 
riod gap and suggests a donor mass around M2 ^0.1 Mq. 
This is surprising. According to standard evolutionary mod- 
els, a CV below the period gap should be faint and have a 
low accretion rate driven primarily by gravitational radia- 
tion (GR). Yet T Pyx's quiescent luminosity and status as 
a RN both imply t hat it ha s a high ac cretion rate of > 10~^ 
M© yr"^ (PattersoireraD[l998 : Selvelli et al]l2008 l). This 
is about two orders of magnitude higher than expected for 
ordinary CVs at this period. 

Assuming that the determination of the photometric 
orbital period is correct, the existence of T Pyx is interest- 
ing for at least two reasons. First, unless we are seeing the 
system in a transient evolutionary state, its lifetime would 
have to be very short r ^ M2/M2 ^ 10 Myrs. This would 
imply the existence of an evolutionary channel leading to 
the fast destruction of at least some short-period CVs. Sec- 
ond, TNRs frequent enough to qualify as RNe are thought 
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to be possible only on high-mass accreting WDs (Mi ^ 1 
M0). Moreover, RNe are the only class of novae in which 
the WD is expected to gain more mass between eruptions 
than it loses during them. This would make T Pyx a strong 
candidate Type la supernova progenitor. 

However, the recent study of the 
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T Pyx was an ordinary cataclysmic variable until it 
erupted as a nova in 1866. This eruption triggered a 
wi nd- driven super soft X-ray pha se (as first suggested 
bvlKnigge. King fc P atterson! |200Q| 1, resulting in an unusu- 
ally high luminosity and accreti on rate. However, unlike in 
the original scenario proposed by Knigg e, King fc PattersonI 
(^00), the supersoft phase is not self-sustaining, so that 
the accretion rate has been declining ever since the 1866 
nova eruption from M ^ 10"'^ Moyr"^ to lO^^Moyr"^ 
As a result, T Pyx has fade d by almost 2 magnitudes 
sinc e the nova eruption (| Schaefer, Pagnotta fc Sharal 
[200s). Based on this, and the fact that T Pyx has al- 
ready passed its mean recurrence time by more than 20 
years, Schaefer, Pagnotta & Shara (2009) argue that T Pyx 
might no longer even be a recurrent nova. If these ideas are 
correct, T Pyx is not a viable SN la progenitor, and its 
remaining lifetime can be substantially longer than a few 
million years. However, if all its ordinary nova eruptions 
are followed by relatively long-lived (> 100 yrs) intervals 
of wind-driven evolution at high M, its secular evolution 
may nevertheless be strongly affected, with significant 
implications for CV evolution more generally (see also 
Section [5]) . 

A key assumption in virtually all of th ese arguments is 
that t he photometric period measured by [Patterson et al.l 
(| 19981 ) is, in fact, the orbital period of the system. So far, 
there has only been one attempt to o btain a spectroscopic 
period for T Pyx, bv lVogt et al.l (|l990l l, who reported a spec- 
troscopic modulation with P — 3.44 hours. Such a long or- 
bital period above the CV period gap would be much more 
consistent with the high accretion rate found in T Pyx. In 
this study, we present the first definitive spectroscopic de- 
termination of the orbital period of T Pyx, showing that 
it is, in fact, consistent with Patterson et al.'s photometric 
period. We also use our time-resolved spectroscopy to esti- 
mate the main system parameters, such as the velocity semi- 
amplitude of the white dwarf (i^l), the mass-ratio (q), the 
masses (Mi and M2) and the orbital inclination (i). Finally, 
we discuss the implications of our results for the evolution 
of T Pyx and related systems. 



2 OBSERVATIONS AND THEIR REDUCTION 

2.1 VLT Multi-fibre Spectroscopy 

Multi-fibre Spectroscopy of T Pyx was obtained during five 
nights in 2004 and 2005 with the GIRAFFE/FLAMES in- 
strument mounted on the Unit Telescope 2 of the VLT at 
ESO Paranal, Chile. The data were taken in the integrated 



field-unit mode. The total field of view in this mode is about 
11.5" X 7.3" a nd thus covers most of T Pyx's 10" diame- 
ter nova shell (Williams| Il982h . We used the fibre system 
ARGUS, which consists of 317 fibres distributed across the 
field, of which 5 are pointing to a calibration unit and 15 
are pointing to sky. The grating order was 4, which gives 
a resolution of R= 12000. The wavelength range was chosen 
between 4501 to 5078 A, so that the emission-line spectrum 
would include the Bowen blend at 4645 - 4650 A and the 
Hell at 4686 A. With this setup, the dispersion is 0.2 A/pix, 
corresponding to about 12.5 kms~^/pix. The full widths at 
half-maximum (FWHMs) of a few spectral lines obtained 
simultaneously with the science spectra from fibres pointing 
to the calibration unit indicate that the spectral resolution 
is about 0.4 A. A log of the observations can be found in 
Tabled 

The initial steps in the data reduction were performed 
using the ESO pipeline for GIRAFFE. The pipeline is based 
on the reduction software BLDRS from the Observatory 
of Geneva. The basic functions of the pipeline are to pro- 
vide master calibration data and dispersion solutions. The 
pipeline also provides an image of the reconstructed field of 
view, which can be used to associate a specific fibre to a 
given object. In order to extract the spectrum from the de- 
sired fibre and to correct for the contribution from the sky 
background and cosmic rays, the output from the pipeline 
was processed further in IRAF. The PSF of our target, T 
Pyx, is spread out over several fibres and 6 single fibre spec- 
tra containing significant target flux were extracted and me- 
dian combined after weighting each spectrum by the mean 
flux in the region 4660 - 4710 A(covering Hell at 4686 A). 
Cosmic rays were removed by first binning the spectra over 
7 pixels (the cosmic rays have a typical width of 2 — 6 pix- 
els). The smoothed spectra were then subtracted from the 
corresponding fibre spectra to only leave the residuals and 
the cosmic rays. Cosmic rays were then removed and the 
residuals added back to the target spectra. Fibres containing 
the sky were extracted and combined to create a master-sky 
spectrum. Finally, the master-sky was subtracted from the 
combined science spectrum. 

2.2 Magellan Long-slit Spectroscopy 

T Pyx was observed again during four nights in March 2008, 
this time on the 6.5 meter Baade telescope Magellan I at 
Las Campanas, Chile (see Table [1] for a log of the observa- 
tions) . Long-slit spectroscopy was obtained with the instru- 
ment IMACS using the Gra-1200-17.45 grating with a 0.9" 
slit, resulting in a dispersion of 0.386 A/pix, corresponding 
to about 25 kms~"^/pix. We estimate that the spectral res- 
olution of these observations is about 1.5 A, as measured 
from the FWHMs of a few spectral lines in the arc-lamp 
spectra. The overall spectra span over four CCDs and cover 
a total wavelength range of 4000 - 4800 A. All frames from 
the four CCDs were treated separately during both the 2D 
and ID reduction steps. 

The spectra were reduced in IRAF using standard 
packages. A master bias produced from combining all bias 
frames obtained during the four nights was subtracted from 
the science frames and the overscan region was used to 
remove the residual bias. A master flat-field, corrected for 
illumination effects was produced. The spectra were then 
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Table 1. Log of the observations (the date is according 
to UT time at the start of the night). 



Date 


Tel/Inst. 


No. of Exp. 


Exp.[s] 


041210 


VLT/GIRAFFE 


28 


180 


041224 


VLT/GIRAFFE 


28 


180 


050127 


VLT/GIRAFFE 


29 


180 


050128 


VLT/GIRAFFE 


28 


180 


050131 


VLT/GIRAFFE 


29 


180 


080316 


MAGELLAN/IMACS 


37 


240 


080317 


MAGELLAN/IMACS 


62 


240 


080318 


MAGELLAN/IMACS 


62 


240 


080319 


MAGELLAN/IMACS 


39 


240 



flat-field corrected and extracted. Line identifications of 
the ThArNe lamp spectra obtained during the nights were 
carried out with the help of NOAO Spectral Atlas Central 
and the NIST Atomic Spectra Database. The spectra were 
cosmic ray rejected using the method described previously 
for the VLT dataset. A flux calibration was done using data 
of the white-dwarf flux-standard star EG274, obtained dur- 
ing the same nights as the target spectra. Also, extinction 
data from the site and flux reference data of EG274 from 
iHamuv, et al.l (|l992h were used. 

Analysis of the reduced VLT and Magellan data were 
carried out using the packages MOLLY and DOPPLER, 
provided by Tom Marsh. 



3 DATA ANALYSIS 

3.1 The Overall Spectrum 

T Pyx has a high-excitation emission-line spectrum that is 
unusual when compared to systems with a low M, but the 
overall spectrum is similar to that for nova-like variables, 
supporting the idea that T Pyx has a high accretion rate. 
Its brightness has been found to be fading for the last cen- 
tury, indicating a decrea se in M, and its 2009 magnitude i n 
blue was B = 15.7 mag (" Schaefer, Pagnotta fc Sharall2009l l. 
Double-peaked lines originating from the disc are seen in the 
spectrum, and the bright disc outshines any spectral signa- 
ture of the primary WD and the donor star. The strongest 
feature is the doubled-peaked Hell line at 4686 A. Double- 
peaked lines identified as Hel at A4713, A4921 and A5015, as 



well as the Balmer lines, are also present (Figure 1(a) and 



Figure 1(b)) 



The Bowen blend at A4640 — A4650 is clearly visible and 
consists of several lines of NIII and CIII. The Bowen blend is 
also seen in other CVs and in low-mass X-ray binaries with 
high-accretion rates. In some cas es, it has been rela ted to 
emission from the donor star feg. ISteeghs fc Casared [2003 1, 
where it is thought to be produced by a fluorescence pro- 
cess as the front side of the donor is strongly irradiated by 
the hot accretor. The Bowen blend was clearly visible in the 
spectrum of T Pyx obtained by Margon k. Deutsch (1998), 
and based on this, we were hoping to detect narrow compo- 
nents associated with the donor. However, no narrow donor 
star features were found in the blend. 

Data analysis was carried out for all the strongest lines 
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Figure 2. Power spectrum of the Hell radial- velocity data, ob- 
tained from the combined VLT and Magellan data sets. The ver- 
tical grey line in the bottom panel represents the best period 
found by Patterson et al. 1998, along with its error. The red 
line represents the recalculated period using the latest photomet- 
ric ephemeris (see Section [3.2. 1|) where the period derivative has 
been taken into account (Patterson et al. in preparation). 



but our final results are based on analysis of the Hell line 
at 4686 A and the Hel line at 4921 A. 



3.2 The Orbital Period 

3.2.1 The Photometric Ephemeris 

The key goal of our study is to obtain a definitive deter- 
mination of T Pyx's orbital period based on time-resolved 
spectroscopy. More specifically, w e wish to test if the 
photometric modulation reported bvl Patterson et al.l (| 19981 ) 
is orbital in nature. The photometric ephemeris is given 
bv lPatterson et al.l (|l998l l as 



Minimum light = 
3.5 X l^-^^E^. 



HJD 2446439.428 + 0.0762233^ + 



The quadratic term is highly significant and indicates 
a very high period derivative of P :^ 9 x 10"^*^ (this corrects 
a typo in the original paper). The period is found to be 
increasing, which is not expected for a system that has not 
yet passed the minimum period. If the photometric signal 
is orbital, the current evolutionary time-scale of the system 
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(a) The normalized and averaged spectra constructed from 140 
individual exposures obtained at the VLT telescope. 
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(b) The flux-calibrated average of 200 spectrum observed with the 
Magellan telescope. 



Figure 1. 



is only r ^ P/P ^ 3 x 10^ yrs. Such a short period change 
time- scale is highly unusual for a CV. 

The Center for Backyard Astronomy (CBA) has 
continued to gather additional timings in the period 1996 
- 2009. This new photometry, along with an updated 
ephemeris, will be published in due course (Patterson et al., 
in preparation). However, since this new data set is much 
more suited to a comparison against our spectroscopic 
observations (since it overlaps in time), Joe Patterson has 
kindly made this latest set of CBA timings available to 
us. Our own prelimary analysis of all the existing timings 
yields a photometric ephemeris of 

Minimum light = HJD(UTC) 2451651.65255(35) + 
0.076227249^(16) +2.546(54) x 10"^^^^ 



The numbers in parenthesis are the uncertainties on 
the two least significant digits, and were estimated via boot- 
strap simulations. Taking the period change into account, 
the average period at the time of the optical spectroscopy 
was P = 0.076228860 =b 3.0 x 10"^ d (corresponding to 
/ = 13.1183912 zb 5.2 x 10"^ c/d). 

This new ephemeris confirms that the period contin- 
ues to increase at a fast rate. More specifically, the current 
estimate of the period derivative is P = (6.68=b0.14) x 10"^*^, 
corresponding to a timescale of P/P ^ 3 x 10^ yrs. It is 
important to take this period derivative into account when 
comparing spectroscopic and photometric periods. 



3.2.2 The Spectroscopic Period 

In order to establish the orbital period spectroscopically, 
we carried out a radial-velocity study usin ^ the double- 
Gaussian method (S chneider fc Youndll980l l. This method 
is most effective in the line wings, which are formed in the 
inner regions of the accretion disc. As discussed further in 
Section [331 we tried a variety of Gaussian FWHM's, as well 
as a range of separations between the two Gaussians. 

The results shown in this section are for a FWHM = 



450 kms~^ and a separation of 300 kms~^, and are based 
on the full combined VLT + Magellan data set, spanning 
the time period December 2004 to March 2008. The radial- 
velocity study was carried out for several spectral lines but 
for the purpose of this section, we will focus exclusively 
on the radial-velocities derived from the strong Hell line 
at 4686 A, since this line p rovides the m ost accurate re- 
sults. The power spectrum (|Scarglelll982i ) of this spectro- 
scopic data set is presented in Figure [2l The top panel, 
covering a broad frequency range, confirms that there is a 
strong signal near the photometrically predicted 13 c/d. No 
signal is found a t P = 3.4 4 hour s, corresponding to the 
period found bv IVogt et al.l (|l990l V The middle and bot- 
tom panels show close-ups of narrower frequency ranges 
around this power spectral peak. Note, in particular, that 
the red line in the bottom panel marks the photometri- 
cally predicted frequency for the epoch corresponding to the 
mid-point of our spectroscopic data set. For comparison, we 
also s how the frequ e ncy of the linear ephemeris determined 
by P atterson eTsU ([l998) for just the 1996-1997 timings, 
along with its error. The spectroscopic and photometric pe- 
riods appear to be mutually consistent, but only if we ac- 
count for the photometry period derivative. More specifically, 
the P- corrected photometric prediction lies extremely close 
to the peak of the most probable spectroscopic alias. 

This apparent consistency can be checked more quan- 
titatively. The orbital frequency predicted by the photom- 
etry for the average epoch of our spectroscopy is fphot — 
13.1183912 ± 5.2 x 10~^ c/d, where the error is once again 
based on bootstrap simulations. Similarly, we have carried 
out a bootstrap error analysis for the frequency determined 
from our spectroscopy. Since we are only interested in the 
agreement between the photometry and the spectroscopy, 
we only considered the spectroscopic alias closest to the pho- 
tometrically determined frequency (this is also the highest 



alias). This yielded fs-^ 



13.118368 ± 1.1 X 10"^ c/d. The 



difference between these determinations formally amounts 
to just under 2a. We consider this to be acceptable agree- 
ment. 

The level of agreement between spectroscopic and pho- 
tometric periods is important. The spectroscopy on its own 
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is sufficient to establish beyond doubt that T Pyx is a short- 
period system with Porh — 1.8 hrs. However, the fact that we 
can estabhsh reasonable consistency with the photometric 
data if (and only if) we account for the photometric period 
derivative suggests that (i) we can trust the much higher 
precision photometric data to provide us with the most ac- 
curate estimate of the orbital period, and (ii) that the large 
period derivative suggested by the photometry is, in fact, 
correct. 



3.3 Trailed Spectrograms and Doppler 
Tomography 

In order to visualize the orbital behavior of the spectral pro- 
files, trailed spectrograms were constructed for Hell, Hel, 
H/3 and the Bowen blend. These spectrograms, shown in Fig- 
ure[3l are phase binned to match the orbital resolution of our 
data (37 bins for the VLT data and 27 bins for the Magellan 
data), and are here plotted over 2 periods. All lines show 
emission features moving from the blue to the red wing, and 
localized emission is seen in the red wing at phase, ~ 0.6. 
Similar structures are seen in the spectrograms for both the 
HeH and the Bowen blend indicating that they originate in 
the same line- forming region, presumably the accretion disc. 
Doppler tomography (Marsh & Home 1988), is an in- 
direct imaging method where the emission-line profiles, de- 
pending on phase, are plotted onto a velocity scale. The 
method was used to visualize, in velocity space, the origin 
of the emission features found in T Pyx. In Figure^ Doppler 
tomograms for the most prominent lines are presented. For 
the reconstruction of the tomograms, the systemic velocity, 
7 was set to zero (see Section [3.4.ip . All lines show asymmet- 
rically distributed emission from the disc. These asymmetric 
features are also seen in the trailed spectrograms (Figure [3]). 
No emission can be connected to the bright spot, the accre- 
tion stream or the donor star. The map of the Bowen blend 
was constructed from a composite of the three lines, NIII at 
A4640.64 and CHI at A4650.1 and A4647.4. 



3.4 A Radial- Velocity Study of T Pyx 

3.4-1 The Systemic Velocity 

We noticed early on in our analysis for both the VLT 
and Magellan data sets that the systemic velocity, 7, ap- 
peared to be varying with time. In order to confirm the 
reality of these variations, we extracted spectra from a 
few VLT/GIRAFFE/FLAMES fibres that were centered on 
knots in the nova shell surrounding T Pyx. The intrinsic 
radial velocities of these large-scale knots are not expected 
to change appreciably over a time-scale of months, so they 
provide a useful check on the stability of our wavelength 
calibration. We found that the shifts seen in the radial ve- 
locity data for the central object are much larger than those 
seen for the knots, suggesting that these shifts cannot be 
explained by instrumental effects alone. However, since we 
have only five nights worth of VLT data spread over two 
months, our data set is too sparse to allow a careful study 
of trends in the apparent 7 velocities. For the purpose of the 
present study, we thus subtracted the mean nightly 7 from 
all radial velocities before carrying out any analysis. This 
should minimize the risk of biasing our results. 




200 400 600 800 1000 1200 1400 

Separation a [km/s] 

Figure 5. Diagnostic diagram for the Hell line, Magellan. None of 
the parameters (7, K and ^o) are converging for any combination 
of the FWHM and separation. 



3.4-2 The Velocity Semi- Amplitude 

In order to establish the velocity semi- amplitude, i^l, of the 
WD, further analysis was pursued on the radial- velocity data 
obtained from the double- Gaussian method. A least squares 
fit was made to the radial velocity curves by keeping the 
orbital period fixed, but allowing the three parameters, the 
systemic velocity, 7, the velocity semi-amplitude, K^ and 
the phase, 00, to vary. Phase zero corresponds to the photo- 
metric phase minimum light. The errors on the input radial 
velocity data, provided in MOLLY, were used to calculate 
the x^ and then rescaled so that y^ — 1. The technique was 
applied to several strong lines in both the VLT and Mag- 
ellan data sets, but the radial velocity data measured from 
the Hell line at 4686 A, and in particular from the Magellan 
dataset, gave the most reliab le fits. The results were plotted 
in a diagnostic d iagram (see Shafter 1983: Thorstensenll 19861 
and lHorne, Wade, fc Szkodv..l98 6V However, as can be seen 
in Figure [5l the key parameters (7, K and 00 ) did not con- 
verge convincingly for any combination of the FWHM and 
separation. Normally, it is thought that the bright spot or 
other asymmetries in the disc are responsible for distorting 
the diagnostic diagram. However, in T Pyx the disc appears 
to be rather symmetric, and no contribution from the bright 
spot is seen (see Figure^]). We conclude that the radial veloc- 
ity curves obtained from the double- Gaussian method are, 
most likely, not tracking the velocity of the primary WD, 
and therefore, i^l, and 00 cannot be estimated reliably us- 
ing this method. 
In order to improve on the results obtained from the double- 
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Figure 3. Trailed spectrograms, phase binned to match the orbital resolution (37 bins for the VLT data and 27 bins for the Magellan 
data), are here plotted over 2 periods. Phases are plotted against velocities (in kms~^). 
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Figure 4. Doppler tomograms for the most prominent lines. The center of symmetry in each map is marked as a dot and the outer disc 
radius is plotted as a circle. Both axis are in kms~^. 



Gaussian method, we used two other techniques to obtain 
an estimate of Kl. First, we measured the velocity center 
of the lines by fitting Gaussians to the individual spectral 
lines. Two Gaussians were fitted to the double-peaked lines, 
keeping the FWHM and the rest wavelength of the line fixed, 
but allowing the peak offsets and the peak strengths to vary. 
Radial velocity curves were then plotted for both blue and 
red wings. As can be seen from the trailed spectrograms, 



Figure [3l the red wing (in particular for the Hell lines) does 
not show the same smooth orbital signal as the blue wing. 
Consequently, only the data from the blue wing was taken 
into account when measuring Kl. Several lines from both 
the VLT and Magellan data sets were analyzed. The data 
were phase-binned to the orbital resolution, velocity-binned 
in a ^ 100 pixel wide region around the line, and the nightly 
mean 7 was subtracted. Hell at A4686 is the strongest line 



The orbital period and system parameters of the recurrent nova T Pyx 7 



and was ultimately used to estimate our best-bet value of 
the velocity semi-amplitude of the WD: Kl — 17.9 ±1.6 
kms~^. Figure [6] shows the radial velocities obtained for the 
Hell line in the Magellan data set. Data presented in black 
are phase-binned to the orbital resolution afforded by our 
time- resolved spectroscopy, while all invidual data points are 
shown in grey. Over-plotted onto the radial velocity curve 
is the best sinusoidal fit, indicating a phasing close to zero, 
00 = —0.03 db 0.03. Figure [7] shows the final trailed spectro- 
grams from the Hell at A4686 in the Magellan data set that 
was used to determine Kl. 

Second, Kl can also be measured from the Doppler to- 
mograms, by locating the center of symmetry of the disk 
emission that dominates the maps. Thi s center should be 
locat ed at V^^ = and Vy = —Kl (eg. ISteeghs fc Casarea 
[2002). After masking out any regions of strong asymmet- 
ric non-disk emission, the center of symmetry was found by 
taking a trial point, smoothing the map azimuthally about 
this point, and then subtracting this smoothed map from 
the real map. If the trial point is far from the center of sym- 
metry, the resulting difference map will contain large asym- 
metric residuals. The optimal center of symmetry is thus the 
trial point for which the difference map exhibits the smallest 
residuals (as measured by the standard deviation of the dif- 
ference map). This procedure yielded estimates for Kl that 
were consistent with those obtained from the Gaussian fits 
(14 < i^l < 19 kms"^). We thus adopt Kl = 17.9 it 1.6 
kms~"^ as our best-bet estimate. 

In closing this section, it should be acknowledged that 
our final estimate of i^l is subject to considerably system- 
atic uncertainties. In particular, it is difficult to rule out 
that whatever is preventing the double Gaussian method 
from converging may also bias the Gaussian fits to the line 
peaks and the center of symmetry of the tomograms. Strictly 
speaking, our estimate of Kl should thus perhaps be viewed 
as a lower limit. However, we are reasonably confident that 
our measurements do trace the orbital motion of the WD 
in T Pyx. This is mainly because the presence of persistent 
double-peaked lines and ring-like structure in the tomograms 
implies that there is a strong accretion disk component to 
the line emission. This disk component must ultimately trace 
the motion of the WD, and our two preferred methods are 
designed to exploit this in as direct a fashion as possible, 
while simultaneously minimizing contamination from asym- 
metrically placed emission regions. 



3.4-3 The Velocity at the Outer Disc Radius 

The peak-to-peak separations of the double peaked emis- 
sion lines can be used to estimate the projected velocity at 
the outer disc radius, v{Rdisk) si^^i- We thus fitted Gaus- 
sians to the phase-binned, double-peaked H/3, Hel and Hell 
lines in the Magellan data set. The peak-to-peak separa- 
tions varies significantly with line but the mean value is 
close to that for Hel. We therefore use the Hel peak-to- 
peak separation for our estimation of v{Rdisk) sini, and 
adopt half the full range as a rough estimate of the asso- 
ciated error, AVpeak-to-peak = 290 lb 26 kms""*^. Smak 
(|l98lh investigated the exact relationship between the true 
v{Rdisk) sini and the measured AVpeak-to-peak, taking into 
account the effects of instrumental broadenings and differ- 
ent disk emissivity distributions. He found that, for a wide 
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Figure 6. Radial velocity curve of the Hell line, Magellan. Data 
in black is phase binned while all data are shown in grey. The 
best sin fit is over-plotted onto the phase-binned data. 
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Figure 8. Doppler tomogram for the Hell line, Magellan. The 
location of the WD, Kl, the outer disc radius as well as the 
location of the donor star, K2, and the accretion stream is plotted 
onto the map. 



range of parameters, v{Rdisk)smi = 0.5 AVpeak-to-peak /u, 
where u = 1.05 =b 0.05. Taking this into account, we find 
v{Rdisk) sini = 138 =b 15 kms~^. 



4 SPECTROSCOPIC SYSTEM PARAMETERS 

If we are willing to assume that our estimates of i^l and 
v{Rdisk) sini are valid, we can use our data to estimate sev- 
eral other system parameters. 
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Figure 7. Trailed spectrograms for the Hell line, Magellan. 



4.1 Mass Ratio 

The method we use to estimate the mass ratio, ^, is similar to 
that described in Warner (1973). The basic idea is that, sub- 
ject to fairly benign assumption, the ratio v{Rdisk) sini/Ki 
should be a function of only q. More specifically, we as- 
sume that the disk is circular, Keplerian, and, since T 
Pyx is a high-M system, that it extends all the way to 
the tidal radius. The d isk radius ca n then be taken as 
Rdisk/a — 0.60/(1 + q) (| Warner! 1 1995l l. These assumptions, 
together with Kepler's third law, yield the expected rela- 
tionship 



v{Rdisk)^ini ^ ^ 2909^^ 
Kl ' a 



(1) 



Using this expression, we estimate a mass ratio of 
^ = 0.20 ±0.03. 



Given a mass ratio and disk radius, it is possible to 
estimate the expected phasing of the bright spot in the 
system, if we make the usual assumption that the bright 
spot lies at the point where the ballistic accretion stream 
fro m the Li point meets the disk edge (see, for exam- 



ple, Patterson, Thorstensen fc Kniggell2008l : [Patterson et al] 
bOOd l . Carrying out this calculation for T Pyx suggests that 
the bright spot could be located at orbital phase cj) ^ —0.03. 
This is interesting, because it is consistent with the negli- 
gible phase offset we have found between photometric min- 
imum light and the red-to-blue crossing of the WD radial 
velocity curve. Thus, based on phasing alone, the bright 
spot could be the source of the photometric modulation. 
However, there are also other interpretations for the photo- 
metric signature. For example, the orbital signal could be 
dominated by a reflection effect, i.e. the changing projected 
area of the irradiated face of the donor as it moves around 
the orbit. This would explain why inferior conjunction of 
the secondary corresponds to minimum light. On the other 



hand, it raises the question why we do not see the signature 
of the irradiated donor in our spectroscopy. As noted above, 
there is no sign of the narrow Bowen fluorescence lines that 
one might expect in the case of strong donor irradiation. 



4.2 Component Masses 

Recurrent novae are expected to harbor a more massive WD 
than most ordinary CVs, and even than most classical novae. 
In a classical nova, the mass of the WD is typically close to 
1 M q , both theoret i cally and o bservationally. Based on this, 
both l Selvelh et all (|2008l ) and ISchaefer. Pagnotta fc Sharal 
((20091) agree on a plausible mass range for the primary 
WD in T Pyx of 1.2 5 - 1.4 M© (using theoretical models 



bv lYaron et al.ll2005l l. From the theoretical point of view, 
a high-mass WD is needed to achieve a sufficiently short 
outburst recurrence timescale. 

In T Pyx, the bright accretion disc outshines any spec- 
tral signature from the donor star, and we can therefore 
not constrain the donor mass spectroscopically. However, 
we can use the p eriod- density re lation for Roche-lobe filling 
secondary stars (lEggleton Il983ll to set some constra i nts. 

As shown by IPatterson, Thorstensen fc Kempl ((20051) 
and iKnigga (|2006l l, the donor stars in ordinary CVs below 
the period gap are inflated by approximately 10% due to 
mass loss, relative to ordinary main sequence stars of the 
same mass. Given T Pyx's peculiar evolutionary state, it is 
not clear if this level of inflation is appropriate for T Pyx, 
and we therefore adopt a conservative range of 0% — 20% 
inflation. In order to estimate the mass of the donor star, 
we thus take the theoretical mai n-sequence mass-rad ius rela- 
tion from the 5 Gyr isochrone of lBaraffe et al.l (|l998l l, adjust 
the stellar radius to account for inflation, and then find the 
secondary mass that yields the correct density for T Pyx's 
orbital period. This yields M2 = 0.14 ± 0.03 M©. Note that 
the inferred donor mass decreases with increasing levels of 
radius inflation. The corresponding mass of the WD is then 
Ml =O.7±O.2M0. 
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Since this estimate of the WD mass is lower than ex- 
pected for a recurrent nova, we can also turn the problem 
around. If the mass of the WD is > 1 M©, as theoreti- 
cal nova models imply (jYaron et al.ll2005l ). the mass of the 
donor becomes > 0.2 M©, if our estimate of the mass ratio 
is correct. 

All of these estimates should, of course, be taken with 
a grain of salt. In particular, they rely on (i) the correct- 
ness of our measured Kl and (ii) the assumption that the 
disk extends all the way to the tidal radius. Despite these 
uncertainties, our calculations highlight an important point: 
it is highly unlikely that the donor star in T Pyx is already 
a brown dwarf. For example, if we retain the assumption of 
maximal disk size, Kl ^ 7 kms~^ would be required in or- 
der for M2 ^ 0.07 M© for any Mi > 1 M©. We thus rule out 
the possibility that T Pyx is a period bouncer in the usual 
sense, i.e. that it had already reached the minimum period 
for ordinary CVs (^70 min theoretically, or :^ 80 min obser- 
vationally), in which case its secondary would now be well 
below the Hydrogen-burning limit [j 



4.3 The Orbital Inclination 

The orbital inclination for T Pyx is thought to be low 
due to the sharp spectral profiles and low radial velocity 
amplitude. IShahbaz et al.l (|l997l l suggested a lower limit 
of the orbital inclination for T Pyx of i ^ 6° based on 
the peak-to-peak separatio n of the Ha line in their spec- 
tra. IPatterson et al.l (|l998l l estimated an inclination of i ^ 
10° — 20° due to the lo w amplitude of the orbital signal, 
and lSelvelh et al.1 (|2008h estimated i - 20° - 30°. 

The inclination, i, can be estimated from our data via 
Kl = vi sini, where vi = {27va)/P x {q/{q + 1)) and a is 
the distance between the two components obtained from Ke- 
pler's HI law. This provides a constraint on the system incli- 
nation of z = 10=b2 degrees, for any reasonable combination 
of the component masses. We stress that the real uncertainty 
on the inclination is bound to be larger, because the formal 
error does not include systematic uncertainties associated 
with, for example, possible bias in our Kl measurement and 
the assumption of tidal truncation in our derivation of q. 



5 SUMMARY AND DISCUSSION 

The main result of our study is the spectroscopic determi- 
nation of T Pyx's orbital period. Port — 1-83 hrs. This con- 
firms that the system is a CV below the period gap and 
implies that its current accretion rate is at least 2 orders 
of magnitude higher than that of an ordinary CV at this 
period. We also find that our spectroscopic orbital period 
is cons istent with the pho tometric ephemeris found for T 
Pyx ([Patterson et al.lll998l , an updated version is given in 
Section 13.2. ip . This means not only that photometric tim- 
ings can be used as a more precise and convenient tracer 
of the orbital motion, but also that the large period deriva- 
tive required by the photometric ephemeris marks a genuine 



^ T Pyx, is, of course, a "period bouncer" in the basic sense that 
its orbital period is currently increasing. 



change in the orbital period of the system. In fact, the spec- 
troscopic data are consistent with the photometric period 
only if the period derivative is accounted for. The period 
derivative obtained here from the latest combined photo- 
metric data (P — 6.7 x 10~ ^ °) is slightly lower than that 
obtained by [Patterson et al.l (|l998l l from data up to 1997 
(P ~ 9 X 10~^°). A de crease in the rate of period c hange 
would be in line with Schaefer, Pagnot ta fc Sharal (|2009|) 
scenario that T Pyx's days as a high-M recurrent nova are 
numbered (at least until its next ordinary nova eruption). 
In any case, the average time-scale for period change found 
across all the photometric data are about 3 x 10^ yrs. 

We have also used our spectroscopic data to obtain es- 
timates of other key system parameters, most notably the 
radial velocity semi-amplitude of the WD, Kl — 17.9 it 1.6 
kms~^, and the mass ratio, q — 0.20 ± 0.03. The latter 
estimate rests on three key assumptions: first, that our 
determination of Kl is correct, second, that our estimate 
'^(^dzs/c) sinz is correct, and, third, that the accretion disk 
around the WD in T Pyx is tidally limited. Taken at face 
value, this relatively high value of the mass ratio implies that 
the donor star in the system is not a brown dwarf. Thus T 
Pyx is not a period bouncer. If we assume that the radius of 
the secondary is 10% di 10% infiated relative to an ordinary 
main sequence star of the same mass, we find that its most 
likely mass is M2 = 0.14 ± 0.03 M©. 

Overall, the physical picture that emerges from 
ou r study is consistent with th e scenario proposed 
bv ISchaefer, Pagnotta fc Sharal (|2009i V In particular, they 
suggest that, prior to the 1866 eruption, T Pyx was an or- 
dinary CV. That eruption then triggered a high-M w ind- 
driven phase, as suggested by Knigge. King fc PattersonI 
((20001) to account for T Pyx's exceptional luminosity. How- 
ever, this phase is not quite self-sustaining, so that T Pyx is 
now fading and perhaps not even a RN anymore. In line with 
this picture, we find that the mass ratio and donor mass we 
derive are not abnormally low for a CV at its orbital period. 
This shows that the present phase of high-M accretion can- 
not have gone on for too long already. The hint of a declining 
period derivative may point in the same direction, but this 
needs to be confirmed. 

Does all this mean that the phase of extraordinarily 
high accretion rates T Pyx is currently experiencing will 
have no lasting impact on its evolution? Not necessarily. 
In a stationary Roche-lobe-filling system, the orbital period 
derivative and total mass-loss rate from the donor (wind loss 
+ mass transfer) are related via 



-Lorb 
-Lorb 



3C - 1 M2 
2 M2 



(2) 



where C is the donor's mass-radius index. In T Pyx, which 
has an increasing orbital period, we will take C — —1/3, 
which is appropriate for adiabatic mass-loss from a fully con- 
vective star (e.g. iKnigge, King fc PattersonI [200Q| 1. We thus 
expect that M2/M2 ^ Porb/Porb, which suggests a typical 
mass- loss rate from the donor of M2 ^^ 5 x 10~^ Mq/jt 
in its current state. (The accretion rate onto the WD can 
be lower than this, since, in the wind-driven scenario, much 
of this mass escapes in the form of an irradiation-driven 
outfiow from the donor.) If every ordinary nova eruption 
in T Pyx is followed by ^ 100 yrs of such high-M2 evo- 
lution, the total mass loss from the donor in the luminous 
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phase is ^ 5 X 10~^ M© between any two such eruptions. 
This needs to be compared to the mass lost from the donor 
during the remaining part of the cycle. If this is driven by 
gravitational radiation, the mas s loss rate from the donor 
will b e M2 :^ 5 X 10~^^ Mo/yr (JKnigge. King k Patterson! 
l200d ). The recurrence time of ordinary nova eruptions fo r 
such a system is on the order of 10^ yrs (|Yaron et al.ll2005l ). 
so the total mass lost from the donor during its normal evo- 
lution (outside the wind-driven phase) is 5 x 10~^ M©. This 
shows that the long-term secular evolution may be dominated 
by its high-M wind-driven phases, even if the duty cycle of 
these phases is very low (e.g. 0.1% for the numbers adopted 
above). 

It is finally tempting to speculate on the relevance of 
"T Pyx- like" evolution for ordinary CVs. At first sight, the 
numbers above suggest that the evolution of a CV caught in 
such a state may be accelerated by about an order of mag- 
nitude. This is interesting, since it could bear on the long- 
standing problem that there are fewer short-period CVs and 
period bou ncers in current samples than theoretically ex- 
pected (e.g. Patterson et al. 1998; Pretorius, Knigge & Kolb 
I2QO7I : IPretorius et al. 2007: Pretorius k Knigge 2008ab). 
Moreover, this channel need not be limited to systems con- 
taining high-mass WDs that are capable of becoming RNe. 
After all, it is not the recurrent nova outbursts that are of 
interest from an evolutionary point of view, but simply the 
existence of a prolonged high-M phase in the aftermath of 
nova eruptions. In CVs containing lower-mass WDs, such 
a phase may still occur, although its evolutionary signifi- 
cance could still depend on the WD mass. For example, the 
interoutburst time-scale is longer for low-mass WDs, and 
the duration of the high-M phase could also scale with WD 
mass. One obvious objection to this idea is that, observa- 
tionally, most nova eruptions are not followed by centuries- 
(or even decades-) long high-M phases. However, this need 
not be a serious issue. Most observed novae are long-period 
systems, so if the triggering of a wind-driven phase requires 
a fully convective donor, most novae would not be expected 
to enter such a phase. It may be relevant in this context 
that at least one other short-period nova - GQ Mus - ex- 
hibited an exc eptionally long post-ou t burst supersoft phase 
of - 10 years (jOegelman et all 1 19931 ). ICreiner et aD (|2003l l 
have also suggested that the duration of the supersoft phase 
in novae may scale inversely with orbital period. 

However, there is another important consequence to the 
idea that the evolution of many short-period CVs is accel- 
erated by "T Pyx-like" high-luminosity phases. The orbital 
period- derivative is positive in the high-M phase, but nega- 
tive during the remaining times of GR-driven evolution. But 
since Porb/Porb — M2/M2 in both phases, the sign of the 
secular (long-term-average) period derivative will generally 
correspond to the phase that dominates the secular evolu- 
tion. So whenever wind-driving dramatically accelerates the 
binary evolution, the direction of period evolution will also 
be reversed. Is this a problem? Perhaps. The recent detec- 
tion of the long- sought period spike in the distribution of CV 
orbital periods (Gansicke et al. 2009) suggests that there is, 
in fact, a reasonably well-defined minimum period for CVs, 
as has long bee n predict ed by the standard model of CV 
evolution (e.g. iKolblll993l V On the other hand, the location 
of the observed spike (:^ 83 mins) is significantly different 
from the expected one (^ 65 — 70 min). Is it possible that 



the observed period minimum corresponds to the onset of 
wind-driving in most CVs? T Pyx, with Porb — 110 min 
would then have to be an outlier, however, perhaps because 
of an unusually high WD mass. 

The idea that T Pyx-like phases may significantly affect 
the evolution of many CVs is, of course, highly speculative, 
and we do not mean to endorse it too strongly. However, 
it highlights the importance of understanding T Pyx: until 
we know what triggered the current high-luminosity state, 
it will remain difficult to assess the broader evolutionary 
significance of this phase. Note that the apparent uniqueness 
of T Pyx is not a strong argument against such significance. 
For example, if the duty cycle of high-luminosity phases is 
^ 0.1%, as suggested by the numbers above, we should not 
expect to catch many CVs in this state. Thus T Pyx could 
be the tip of the proverbial iceberg. 
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